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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  under  Program  Element  65807F. 
The  results  of  this  research  were  obtained  by  ARO,  Inc.,  AEDC  Division  (a  Sverdrup 
Corporation  Company),  operating  contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force 
Station,  Tennessee,  under  ARO  Project  No.  P32A-K9A.  A.  F.  Money  was  the  Air  Force 
project  manager.  Analysis  of  the  data  was  completed  on  September  30,  1977,  and  the 
manuscript  was  submitted  for  publication  on  January  11,  1978. 
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1.0  INTRODUCTION 


The  simulation  of  aircraft  motion  through  analytical  techniques  is  becoming  a very 
important  tool  in  the  development,  testing,  and  operational  phases  of  modern  fighter 
aircraft  programs.  A pilot  in  the  loop  fixed-base  simulators,  which  in  the  past  have  been 
used  primarily  for  pilot  training  and  proficiency  checks,  are  now  being  applied  to  the 
development  and  testing  phases  of  new  fighter  aircraft.  An  example  of  this  is  the  NASA 
differentially  maneuvering  simulator  which  has  recently  been  used  to  study  the 
effectiveness  of  an  automatic  control  system  on  an  advanced  fighter  aircraft  at  high 
angles  of  attack  (Ref.  1).  By  using  the  simulator,  a more  economical  and  safer 
development  and  evaluation  program  was  conducted  than  was  possible  through  flight 
testing. 

In  the  low  angle-of-attack  unstalled  flight  regime,  the  confidence  level  in  motion 
simulation  is  high.  As  the  aircraft  angle  of  attack  increases  to  the  stall/spin  flight  regime, 
the  confidence  level  diminishes  correspondingly.  This  degradation  in  confidence,  resulting 
from  poor  "before-the-fact"  simulation,  is  not  caused  by  an  inadequacy  in  the  aircraft 
equations  of  motion  in  stall/spin  flight  but  results  from  improper  or  inadequate  definition 
of  the  aircraft  aerodynamics  in  this  regime. 

New  experimental  test  apparatus  in  the  form  of  dynamic  balances  (Ref.  2)  and 
rotational-balance  hardware  at  NASA  Langley  and  Ames  are  being  developed  to  better 
define  aircraft  aerodynamics  in  stalled  and  spinning  flight.  The  rotational  balance 
apparatus  measures  the  aerodynamics  of  an  aircraft  under  steady  rotating  conditions.  This 
accounts  for  the  nonlinear  effect  of  rotation  rate.  Preliminary  analysis  indicates  these 
data  are  necessary  for  before-the-fact  simulation  of  an  aircraft  spinning  motion.  Little 
work  has  been  done  to  define  what  dynamic  oscillatory  components  should  accompany 
these  rotational-balance  data  in  defining  data  matrices  for  spinning  flight.  Generally,  the 
oscillatory  direct  and  cross-damping  derivatives  are  included  in  spin  simulation  because  of 
their  proven  importance  in  the  nonspinning  regime.  The  importance  of  dynamic 
cross-coupling  and  acceleration  derivatives  in  spinning  flight  is  not  clearly  defined. 

The  Arnold  Engineering  Development  Center  (AEDC)  is  presently  conducting 
research  to  define  what  aerodynamic  parameters  are  necessary  for  achieving  good  before- 
the-fact  aircraft  motion  simulation  in  the  stall/spin  flight  regime.  It  is  necessary  to 
define  these  parameters  such  that  future  experimental  test  apparatus  may  be  designed  for 
their  measurements. 


The  subject  analysis  defines  the  importance  of  aircraft  cross-coupling  and 
acceleration  damping  derivatives  in  the  analytical  simulation  of  aircraft  spinning  motion. 
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Derivatives  investigated  in  the  analysis  are  Cmp.  Cmr,  Cfiq.  C„q,  C„a,  and  Cgj.  Two  types 
of  aircraft  are  utilized  in  the  analysis.  Changes  in  one  of  the  predominate  spin  modes  of 
each  aircraft  are  investigated  with  both  individual  and  simultaneous  variations  of 
cross-coupling  and  acceleration  derivatives.  Time  histories  of  the  aircraft  spinning  motion 
are  provided  using  a six-degrce-of-freedom  (6-DOF)  nonlinear  computer  program.  The 
analysis  addresses  only  the  developed  spinning  mode  of  flight:  attempts  to  analyze  the 
incipient  or  recovery  phases  of  the  spin  are  not  included. 

2.0  BACKGROUND 

Prior  to  the  1960's,  aircraft  dynamic  cross-coupling  and  acceleration  derivatives  were 
generally  considered  insignificant  in  fighter  aircraft  motion  simulation  and  dynamic 
stability  analysis.  These  assumptions  were  good,  primarily  because  aircraft  of  that  era 
were  operating  at  moderate  angles  of  attack  where  both  cross-coupling  and  acceleration 
derivatives  possessed  small  values.  In  the  mid  1 960’s,  it  became  necessary  for  these  same 
aircraft  to  extend  their  operating  envelopes  to  high  angles  of  attack  under  combat 
conditions.  As  a result  of  this  high-risk  operating  environment,  aircraft  were  pushed 
beyond  stalled  and  spinning  flight.  Attempts  to  analytically  simulate  the  poststall 
gyration  and  spinning  flight  experienced  by  the  aircraft  were  generally  unsuccessful.  When 
a successful  correlation  was  achieved,  it  was  normally  an  after-the-fact  correlation 
resulting  from  many  attempts  of  parameter  variations  to  achieve  a suitable  combination 
of  aerodynamics.  It  became  apparant  that  a better  representation  of  the  aerodynamics  of 
an  aircraft  would  be  necessary  for  achieving  a before-the-fact  analytical  simulation  of 
aircraft  stalled  and  spinning  motion. 

An  attempt  to  isolate  some  of  the  important  aerodynamic  parameters  in  aircraft 
spinning  flight  was  undertaken  by  Bihrle  in  1970  (Ref.  3).  Dynamic  cross  and 
acceleration  derivatives  were  evaluated  against  spin  sensitivity  in  the  analysis.  Results 
indicated  the  spinning  motion  to  be  insensitive  to  the  acceleration  derivatives  and  often 
sensitive  to  cross  derivatives.  Since  Bihrle 's  work,  it  has  been  shown  experimentally  by 
NASA  (Ref.  4)  and  Orlik-Ruckemann  of  the  NAE  in  Canada  (Ref.  5)  that  the  magnitude 
of  the  acceleration  and  cross  derivatives  of  aircraft  at  high  angles  of  attack  is  much 
greater  than  Bihrle  had  assumed.  Also,  the  sign  of  these  derivatives  often  changes  at  the 
higher  angles.  In  addition,  Ref.  5 shows  the  aircraft  cross-coupling  derivatives  to  have 
values  equal  in  magnitude  to  direct  derivatives. 

3.0  METHOD  OF  ANALYSIS 

The  analytical  tool  used  in  this  dynamic  sensitivity  analysis  was  a 6-DOF  nonlinear 
computer  program.  The  program  was  originally  formulated  by  North  American  Rockwell 
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(Ref.  6)  using  a fourth-order  Runge-Kutta  integration  algorithm  with  a fixed  integration 
size.  The  program  input  and  aerodynamic  modules  have  since  been  modified  for 
adaptation  to  this  study.  The  program  was  modified  to  run  interactively  with  a graphics 
terminal  to  expedite  the  analysis.  The  6-DOF  equations  of  motion  describing  the  aircraft 
motion  are  rigid-body  equations  referenced  to  a body-fixed  axis  system  shown  in  Fig.  1. 
The  basic  equations  are  as  follows: 


Forces 


Moments 


u = 

rv  - 

qw  — 

q sin  6 + * 

lx 

IN 

IN 

• 

V = 

pw  - 

ru  -r 

q cos  0 sin  $ + 

in 

W - 

qu  - 

pv  + 

q cos  0 cos  <f>  4 

. Iz  + Iz 

m m 


p - 


k - u 

_i L or 


qr  I X Z (r  + p q)  + _X 


<1  - 


1 IX  pr  )■  ixji  (r2  - p2)  + _V  + + _JTC 


U - 1, 


M v M ., 


>Y 


'y  ]y 


f = !l  pq  ixjl  (p  - qr)  + ^3  + 

1 z 1 z 1 z 1 z 


The  external  forces  and  moments  (Fx,  Fy , Fz,  Mx , My,  Mz)  in  the  equations  are 
comprised  of  aerodynamic  coefficients  representative  of  the  aircraft.  Formulation  of  the 
aerodynamic  mathematical  models  used  in  representing  each  of  these  aircraft  is  presented 
in  Appendix  A.  It  should  be  noted  that  the  aerodynamic  data  matrix  for  the  attack 
aircraft  consists  of  static  and  dynamic  oscillatory  data,  whereas  the  data  matrices  for  the 
fighter  bomber  have  static,  dynamic  oscillatory,  and  rotational-balance  data.  All 
aerodynamic  data  used  in  the  mathematical  model  of  Appendix  A are  experimental  and 
are  discussed  in  more  detail  in  Section  4.2.  The  method  used  in  implementing  the 
rotational  balance  data  for  the  fighter  bomber  is  outlined  in  Appendix  B. 

Additional  basic  equations  used  in  the  subject  analysis  are  given  below: 


tan  ^ H 

IV 

j8  = sin  1 

N OO  ' 

= yj  u2  4 v2  + iv2 

to 

ii 

< 

2 ° 

- (|  + r 

up  i vq  wr 


7 


AEDC-TR-78-12 


An  explanation  of  the  steady-state  rotational  rate  vector  (£2SS)  is  given  in  Appendix  B. 
Time  histories  of  the  aircraft  spinning  motion  were  generated  using  this  6-DOF  computer 
program.  Angular  motion  changes  about  the  aircraft  center  of  gravity  and  changes  in  the 
center-of-gravity  path  with  either  individual  or  simultaneous  variations  of  the  derivatives 
were  used  to  ascertain  the  significance  of  the  derivatives  in  spinning  motion  simulation. 

4.0  TECHNICAL  DATA 

4.1  AIRCRAFT  CONFIGURATIONS 

Two  aircraft  configurations  were  chosen  for  the  spinning  motion  sensitivity  study. 
The  selection  of  these  aircraft  was  influenced  by  previous  dynamic  sensitivity  studies 
(Ref.  8)  conducted  at  AEDC  using  the  same  configurations.  Also,  it  was  felt  that  aircraft 
should  be  used  for  which  a general  knowledge  of  the  spin  characteristics  exists.  Aircraft 
similar  to  each  of  the  selected  configurations  have  been  involved  in  stall/spin  test 
programs  (Refs.  9 and  10)  identifying  their  basic  spin  characteristics. 

The  attack-type  configuration,  a single  engine  swept  wing  aircraft,  is  shown  in  Fig.  2 
with  overall  basic  dimensions.  The  fighter  bomber  type,  a twin  engine  and  swept  wing 
aircraft,  is  shown  in  Fig.  3,  along  with  overall  basic  dimensions. 

4.2  BASIC  AERODYNAMIC  AND  INERTIA  DATA 

The  aerodynamic  data  matrices  used  in  modeling  the  attack  aircraft  were  obtained 
from  the  Navy  in  an  unpublished  form.  The  fighter  bomber  data  matrices  were 
formulated  from  data  presented  in  Refs.  1 1 and  1 2.  The  basic  matrices  are  a combination 
of  static  and  dynamic  oscillatory  aerodynamic  data.  Data  are  input  to  the  matrices  in 
table  look-up  form  as  functions  of  the  variables  shown  in  the  equations  of  Appendix  A. 
The  fighter  bomber  data  matrix  also  includes  rotational  balance  data,  previously  used  in 
Ref.  7,  as  an  increment  to  the  basic  aerodynamics.  The  attack  aircraft  data  matrix  does 
not  include  rotational  data.  Neither  of  the  data  matrices  has  been  corrected  for  Reynolds 
number,  flexibility,  or  inertia  effects  experienced  in  full-scale  flight. 

Mass  and  inertia  characteristics  of  the  aircraft  configurations  aTe  presented  in  Table 
1.  Although  the  magnitudes  of  the  mass  and  inertia  of  the  two  aircraft  are  considerably 
different,  their  mass  distributions  along  the  reference  axis  are  similar.  Figure  4 presents 
mass  distributions  of  several  modem  fighter  aircraft  in  the  military  inventory.  With  the 
exception  of  the  F-5,  all  of  these  aircraft  have  similar  mass  distributions  along  the 
reference  axis  similar  to  the  aircraft  used  in  the  project  investigation. 
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5.0  RESULTS  AND  DISCUSSION 

5.1  BASELINE  SPINNING  MOTION 

The  approach  used  for  conducting  the  dynamic  sensitivity  study  was  to  first 
establish  an  equilibrium  baseline  spinning  motion  with  which  other  spinning  motions 
could  be  compared  as  derivative  variations  were  made.  The  baseline  spins  used  for  the 
aircraft  are  shown  in  Figs.  5 and  6. 

The  simulated  spinning  motion  of  the  attack  aircraft,  Fig.  5,  is  classified  as  a steep 
smooth  spin.  The  motion  was  initiated  from  trimmed  1-g  flight  with  control  inputs  of 
stick  back  (5K  = -25  deg),  ailerons  neutral  (5A  = 0 deg),  and  rudder  with  (5k  = -25).  The 
spin  was  entered  at  an  altitude  of  21,600  ft.  Engine  thrust  remained  constant  throughout 
the  spin  at  6,400  lb. 

The  simulated  spinning  motion  of  the  fighter  bomber  aircraft,  Fig.  6,  is  classified  as 
Hat  and  compares  closely  with  flat  spins  experienced  in  flight  tests  (Ref.  10).  The  motion 
was  initiated  by  inputting  conditions  (angle  of  attack,  sideslip,  rotational  rates,  etc.)  at 
time  zero  equal  to  the  equilibrium  flat  spin  parameters  of  flight  data  in  Ref.  10.  By 
initiating  the  spin  in  this  manner,  prospin  controls  were  not  required  for  entering  or 
maintaining  the  spin.  Lateral  directional  control  surfaces  remained  in  the  zero  position 
and  the  elevator  remained  at  3 deg  (5^  = 3 deg).  The  spin  was  initiated  at  an  altitude  of 
12,500  ft.  Engine  thrust  was  maintained  near  zero  throughout  the  maneuver. 

The  dynamic  derivatives  (C^p,  Clllr.  Cfcq,  Cnq,  Cnp-.  and  Qy)  investigated  in  the 

dynamic  sensitivity  study  were  each  given  nominal  values  in  the  above  baseline  spins.  The 
values  used  were  selected  from  experimental  data  (Refs.  4 and  5)  obtained  on 
airplane-like  configurations  and  are  presented  in  Table  2.  The  data  presented  are  for  angle 
of  attack  of  approximately  37  deg  which  corresponds  closely  to  equilibrium  spin  angle  of 
attack  of  the  attack  aircraft.  The  equilibrium  angle  of  attack  ol  the  baseline  spin  for  the 
fighter  bomber  aircraft  is  near  75  deg.  Because  no  experimental  data  exist  for 
cross-coupling  derivatives  in  this  angle-of-altack  range,  it  was  decided  to  use  the  37-deg 
data  of  Table  2 for  the  baseline  spins  for  both  types  of  aircraft. 

5.2  DYNAMIC  DERIVATIVE  VARIATIONS 

The  dynamic  derivatives  under  investigation  were  input  to  the  baseline  spinning 
maneuvers  as  constant  values  at  a time  when  the  aircraft  spinning  motion  had  reached 
equilibrium.  For  the  attack-type  aircraft,  the  spin  equilibrium  was  defined  to  occur  at  20 
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sec  into  the  maneuver,  the  corresponding  time  for  the  fighter  bomber  type  was  defined 
to  be  10  sec.  In  most  cases,  the  derivatives  were  input  individually  at  these  times  and 
remained  constant  throughout  the  remaining  motion.  This  allowed  the  derivative  to  be 
isolated  in  a manner  such  that  its  importance  in  the  aircraft  spinning  motion  could  be 
readily  evaluated.  For  a limited  number  of  cases,  several  of  the  derivatives  were  input 
simultaneously  at  the  above  times  to  provide  a more  realistic  picture  of  what  is  occurring 
in  flight. 

Experimental  data  of  Refs.  4 and  5 provided  a guideline  for  selecting  the  maximum  and 
minimum  values  of  the  derivatives  to  be  input  in  the  baseline  maneuvers.  The  values  selected 
for  each  derivative  are  presented  in  Table  3.  It  should  be  noted  that  the  range  of  experimen- 
tal values  used  in  Table  3 was'  not  obtained  in  a spinning  flight  environment  with  associated 
higher  angles  of  attack  (above  40  deg)  and  the  presence  of  curvilinear  flow.  Also,  the  cross- 
coupling derivative  data  obtained  from  Ref.  5 were  measured  on  a cone-wing  body,  not  an 
actual  airplane  configuration. 

Cm  p Derivative 

Little  is  known  concerning  the  magnitude  of  the  cross-coupling  derivative  Cmp.  The 
experimental  work  performed  by  Orlik-Ruckemann  presented  in  Ref.  5 addresses  only  the 
cross-coupling  derivatives  associated  with  pitching  and  yawing  of  the  aircraft  and  does  not 
include  Cm  resulting  from  a rolling  motion.  Because  these  works  (Ref.  5)  represent  the 
only  known  published  literature  of  measured  aircraft  cross-coupling  derivatives,  the 
relative  magnitude  of  Cmp  continues  to  remain  in  question.  Therefore,  in  the  present 
sensitivity  study  it  was  elected  to  use  Cmp  maximum  and  minimum  values  identical  to 
those  selected  from  Ref.  5 for  the  derivative  Cmr,  -2  to  2 per  radian. 

The  time  histories  for  the  Cm  p variations  compared  to  the  baseline  spinning  motion 
in  which  Cmp  possessed  a nominal  value  (Table  2)  of  zero  is  presented  in  Fig.  7.  For  both 
aircraft,  the  negative  Cm  p coupled  with  a large  negative  roll  rate  (p)  produces  a driving 
pitching  moment  resulting  in  large  longitudinal  motion  changes  with  accompanying 
smaller  latcral/directional  changes.  For  the  fighter  bomber  aircraft,  the  increased 
oscillations  in  the  longitudinal  plane  of  motion  eventually  result  in  the  aircraft  pitching 
out  of  its  equilibrium  flat  spin  and  entering  a form  of  poststall  gyration  at  the  time  of 
approximately  37  sec.  Note  that  in  Fig.  7 the  mean  value  of  the  roll  and  yaw  rates  (p 
and  r)  are  large  compared  to  the  pitch  rate  (q),  which  is  characteristic  of  a spinning 
aircraft.  When  the  large  p lateral  rate  is  coupled  to  even  a small  cross-coupling  derivative 
such  as  Cm  , the  resulting  aerodynamic  pitching  moment  can  have  a significant  effect  on 
the  aircraft  spinning  motion. 
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Cmr  Derivative 

The  second  cross-coupling  derivative  directly  affecting  the  aircraft  longitudinal 
motion  is  the  Cm  f derivative.  As  in  the  case  of  the  Cm  p derivative,  the  pitching  moment 
resulting  from  Cmi  may  be  caused  by  the  high  yaw  rates  r associated  with  spinning 
motion.  Time  histories  generated  with  Cmr  values  of  -2  and  2 per  radian  are  shown  along 
with  the  standard  baseline  motion  in  Fig.  8.  For  the  attack  aircraft,  the  motion  trends 
resulting  from  the  Clllr  variation  are  very  similar  to  these  for  the  CItlp  variation  in  Fig.  7. 

These  similar  trends  are  expected  since  the  ranges  of  both  derivatives  are  identical  and 
the  p and  r rates  arc  of  the  same  sign  and  close  in  magnitude. 

The  trends  for  the  fighter  bomber  spinning  motion  shown  in  Figs.  7 and  8 would  be 
expected  to  be  similar  with  the  same  derivative  variations,  but  they  obviously  disagree. 
The  tighter  bomber  baseline  motion  is  insensitive  to  the  Cm  r derivative  variation.  This 
insensitivity  is  a result  of  the  method  by  which  the  motion  simulation  was  generated.  As 
pointed  out  in  Appendix  A.  the  fighter  bomber  aerodynamic  data  matrix  incorporates 
rotational  balance  data  in  the  manner  described  in  Appendix  B and  Ref.  7.  To 
incorporate  the  rotational  data,  the  total  p,  q,  and  r rates  were  divided  into  oscillatory 
and  steady-state  components.  The  oscillatory  and  steady-state  components  are  in  turn 
applied  to  their  corresponding  derivatives  in  the  equations  of  motion.  Since  the  baseline 
spin  for  the  Fighter  bomber  is  flat  (a  ^ 75  deg),  the  steady-state  component  of  the  yaw 
rate  is  approximately  equal  to  r (see  Appendix  B).  The  oscillatory  component  of  r is  in 
turn  near  zero,  meaning  that  the  resulting  pitching-moment  contribution  from  Cmr  is 
near  zero,  thereby  resulting  in  the  insensitivity  of  the  motion  to  Cmr. 

When  using  rotational  balance  data  (as  in  Appendix  B)  for  spin  simulation,  all 
dynamic  oscillatory  derivatives  caused  by  yaw  rates  become  increasingly  insignificant  as 
angle  of  attack  increases.  For  conventional  spin  simulation  (no  rotational  data),  the  Cmr 
derivative  is  significant  because  of  high  yaw  rates  experienced  in  the  spin. 

Cgq  Derivative 

The  rolling  moment  of  an  aircraft  resulting  from  pitching  in  the  longitudinal  plane  is 
described  by  the  cross-coupling  derivative  Cgq.  Time  histories  of  the  aircraft  spinning 
motion  with  variations  of  this  derivative  are  given  in  Fig.  9.  Although  the  magnitude 
of  q in  the  baseline  spin  is  small  for  the  attack  aircraft,  the  Cfiq  and  q combination 
produces  significant  effects  on  the  aircraft  spinning  motion.  The  negative  C£q  actually 
forces  the  aircraft  out  of  the  spin  indicating  the  aircraft  motion  is  very  sensitive  to  this 
derivative  when  considering  the  magnitude  of  the  q rate.  This  conclusion  is  further 
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verified  when  viewing  the  results  for  the  fighter  bomber  aircraft.  Here,  the  baseline  q rate 
is  oscillatory  and  of  larger  magnitude  than  observed  in  the  attack  aircraft  data.  The  result 
is  a disruption  of  the  fighter  bomber  spin  almost  immediately  after  the  derivative  is 
injected  (t  = 10  sec)  to  the  baseline  motion.  The  data  shown  in  Fig.  9 document  that  the 
Cgq  derivative  is  the  most  important  of  the  cross-coupling  derivatives  in  the  aircraft 
spinning  motion  investigated. 

Cnq  Derivative 

The  second  cross-coupling  derivative  resulting  from  the  pitching  moment  of  an 
aircraft  is  the  yawing  moment  caused  by  pitch  rate  Cnq.  The  motion  sensitivity  to  the 
0,  derivative  is  presented  in  Fig.  10  for  the  two  aircraft.  It  should  be  noted  that  the 
maximum  C„q  value  of  2 per  radian  corresponds  to  the  nominal  value  used  in  the 
baseline  maneuver;  therefore,  the  two  maneuvers  are  identical.  This  results  from  the  fact 
that  the  angle  of  attack  at  which  the  nominal  values  of  Table  2 were  selected  is  the  same 
angle  at  which  C„q  experiences  its  maximum  value.  The  minimum  value  of  C„q  (-2  per 
radian)  has  a small  effect  on  the  motion  of  both  aircraft  when  compared  to  the  baseline. 
The  primary  motion  difference  is  a small  increase  in  the  spin  rates  of  both  aircraft  as 
seen  in  the  r values.  The  small  changes  in  frequency  of  the  motion  for  the  fighter  bomber 
are  considered  to  be  insignificant.  The  C„q  derivative  appears  to  be  the  least  significant 
of  all  the  cross-coupling  derivatives  in  the  aircraft  spinning  motion  investigated. 

Cgp*.  Cnp-  Derivative 

In  addition  to  the  cross-coupling  derivatives,  the  lateral/directional  acceleration 
derivatives  Cfy  and  Cn^  have  also  been  examined.  Recent  work  presented  in  Ref.  4 
indicates  that  these  derivatives  at  high  angles  of  attack  (angles  in  the  spin  regime)  may 
acquire  values  equal  in  magnitude  to  the  lateral/directional  rate  derivatives.  The 
significance  of  the  acceleration  derivatives  of  this  magnitude  on  aircraft  spinning  motion 
is  shown  in  Figs.  11  and  12  for  Cn^  and  C^,  respectively.  The  values  for  the  derivatives 

are  representative  of  those  obtained  experimentally  in  Ref.  4.  The  spinning  motion  for 
the  attack-type  aircraft  shows  significant  effects  caused  by  the  (1  variations.  The  smooth 
baseline  spin  becomes  oscillatory  with  the  input  of  the  negative  Cn(3\  Likewise,  as  shown 
in  Fig.  12,  the  effect  of  the  positive  Cg£  causes  the  aircraft  to  pitch  out  of  the  spin. 

In  both  Figs.  11  and  12  the  acceleration  derivatives  are  seen  to  have  a small 
influence  on  the  motion  of  the  fighter  bomber  aircraft.  Small  changes  in  frequency  are 
noted  but  the  overall  nature  of  the  spin  remains  the  same.  It  should  be  noted  that  the 
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manner  in  which  the  rotational  aerodynamic  data  are  included  in  the  fighter  bomber 
simulation  should  not  affect  the  influence  of  the  j3  derivatives  as  was  the  case  for  the 
aforementioned  Cnil  derivative.  The  acceleration  derivative  variation  indicates  (1)  the 
acceleration  derivatives  can  be  important  in  aircraft  spin  simulation  for  particular  aircraft 
and  (2)  the  acceleration  derivatives  are  configuration  dependent,  and  general  statements 
as  to  their  importance  may  be  misleading. 

Simultaneous  Variations 

The  major  portion  of  the  subject  sensitivity  study  was  conducted  using  individual 
dynamic  derivative  changes  with  all  remaining  derivatives  held  constant.  This  method,  in 
effect,  isolates  the  derivative  being  investigated  such  that  motion  changes  with  different 
derivative  values  can  be  ascertained.  In  actual  flight,  it  is  impossible  to  isolate  a single 
derivative  in  this  manner;  the  magnitudes  of  all  derivatives  normally  vary  simultaneously 
as  the  derivative  in  question  is  changed.  An  attempt  to  more  accurately  simulate  this 
simultaneous  variation  is  accomplished  by  First  selecting  two  sets  of  dynamic  data  as  a 
function  of  the  independent  aerodynamic  variable  angle  of  attack.  These  data  are 
presented  in  Table  4 for  angles  of  attack  of  1 7 and  26  deg.  The  angles  and  corresponding 
derivatives  selected  are  those  from  Ref.  5 for  Cjq  maximum  and  minimum  values  of  2 
and  -2  per  radian.  The  data  in  Table  4 are  in  turn  input  into  the  aircraft  spinning  motion 
together,  rather  than  individually.  The  time  at  which  the  data  are  input  remains  the  same 
as  in  the  individual  cases. 

Time  histories  presented  in  Fig.  13  reflect  the  effect  of  simultaneous  derivative 
changes  on  the  aircraft  maneuvers.  Both  aircraft  experience  significant  changes  in  then- 
spinning  motion  with  input  of  the  derivatives.  For  the  attack  aircraft,  the  motion  trends 
are  similar  to  those  in  Fig.  9 where  Cgq  was  varied  individually.  In  this  case,  the  aircraft 
motion  is  predominated  by  the  Cj>q  derivative  even  though  the  other  derivatives  also 
changed.  For  the  fighter  bomber  aircraft,  the  simultaneous  variations  caused  the  aircraft 
to  pitch  out  of  the  spin  and  enter  a recovery  dive.  These  motions  are  representative  of 
those  in  Fig.  9 when  Cgqis  -2  per  radian  but  not  for  the  2 per  radian  data.  Undoubtedly, 
the  smaller  derivatives  (Table  4)  predominated  over  the  Cgq  value  of  2 per  radian  in  the 
fighter  bomber  motion.  It  is  evident  from  the  results  of  the  fighter  bomber  that  the 
smaller  derivatives  can  have  a significant  effect  on  the  aircraft  motion  when  acting 
together.  Also,  conclusions  based  on  individual  derivative  variations  can  be  significantly 
different,  depending  on  the  values  of  the  baseline  derivatives. 

This  study  has  shown  the  need  for  dynamic  cross-coupling  and  acceleration  derivatives 
in  simulation  of  aircraft  spinning  motion  when  the  derivative  magnitudes  reach  those  given 
in  Table  3.  The  final  determination  of  the  importance  of  these  derivatives  will  surface  when 
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their  true  magnitudes  are  documented  experimentally  for  an  aircraft  configuration  in  a 
realistic  flight  regime.  If,  under  these  conditions,  they  possess  values  approaching  those  in 
Table  3,  then  aerodynamicists  should  consider  including  them  in  spinning  motion  simula- 
tion. The  design  of  future  dynamic  test  apparatus  should  include  the  capability  of  measuring 
these  derivatives  such  that  their  true  magnitudes  can  be  documented  for  an  aircraft  con- 
figuration at  high  angles  of  attack. 

6.0  CONCLUDING  REMARKS 

As  a result  of  this  analysis  for  spinning  aircraft,  the  following  observations  and  conclu- 
sions are  offered  based  on  the  assumption  that  the  derivatives  investigated  can  acquire 
values  of  the  magnitude  given  in  Table  3. 

1.  Because  of  large  r and  q rotational  rates  associated  with  spinning  flight, 
the  pitching  moments  generated  by  cross-coupling  derivatives  Cmt  and 
On  p may  be  equal  or  greater  than  those  generated  by  the  primary 
derivative  Cm  and  q rate  combination. 

2.  All  cross-coupling  and  acceleration  derivatives  Cmp,  Cmr,  Cgq,  Cnq,  Cn(J-, 
and  Cjy  can  produce  significant  effects  on  the  aircraft  spinning  motion. 

3.  Spinning  motion  sensitivity  to  the  dynamic  cross-coupling  and  acceleration 
derivatives  is  configuration  dependent. 

4.  The  sensitivity  of  an  aircraft's  spinning  motion  to  an  individual  dynamic 
derivative  variation  may  vary  as  a function  of  the  baseline  derivatives  used 
in  the  investigation. 

5.  The  final  determination  of  the  importance  of  the  cross-coupling  and 
acceleration  derivatives  will  remain  in  question  until  their  magnitudes  are 
experimentally  measured  for  an  aircraft  configuration  in  a realistic  flight 
regime;  experimental  test  apparatus  should  be  developed  for  ascertaining 
these  magnitudes. 
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Figure  1.  Aircraft  body-axis  system. 
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Figure  4.  Aircraft  mass  distribution. 
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Figure  7.  Concluded. 
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Figure  10.  Concluded. 
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° Attack  Baseline  Spin 
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a.  Attack-type  aircraft 
Figure  12.  Spin  sensitivity  to  variation. 
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b.  Fighter  bomber-type  aircraft 
Figure  12.  Concluded. 
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Table  1.  Aircraft  Physical  and  Mass  Characteristics 

Attack-Type  Aircraft 


Mass  719  slugs 

Ix  15,927  slug-ft2 

IY  64,792  slug-ft2 

Iz  75,976  slug-ft2 

lxz  3,885  slug-ft2 

S 375  ft2 

b 38.73  ft 

c 10.84  ft 

eg  30  percent  MAC 


Fighter  Bomber-Type  Aircraft 

1147.5  slugs 
22,600  slug-ft2 
163,100  slug-ft2 

182.000  slug-ft2 

2 

5,450  slug-ft 
538.3  ft2 

38.4  ft 

16.04  ft 

33.0  percent  MAC 


Table  2.  Cross-Coupling  and  Acceleration  Derivatives 
for  Baseline  Spin  (Body  Axis) 

c . =0.25 

ne 

C„.  = -0.35 

C = -0.75 
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c = 2.0 
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q 


37 


AEDC-TR-78-12 


Table  3.  Dynamic  Derivative  Values  for  Individual 
Variations  (Body  Axis) 


Max 

Min 

C . 

0 

1.0 

-0.5 

cr 

1.0 

-1.0 

C 

m 

r 

2.0 

-2.0 

C 

m 

P 

2.0 

-2.0 

C* 

q 

2.0 

-2.0 

c 

n 

q 

2.0 

-2.0 

Table  4.  Dynamic  Derivatives  Used  in  Simultaneous 
Variations  (Body  Axis) 

Case  for  Angle  of  Attack  =17  deg  Case  for  Angle  of  Attack  = 26  deg 
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c . 
n6 

= 

0.71 

= 

-0.5 

cr 

3 

= 

-0.95 

C 

m 

r 

0.4 

C 

in 

r 

ss 

-0.2 

C 

m 

= 

0 

C 

m 

= 

0 

P 

P 

-2.0 

=5 

2.0 

q 

q 

C 

n 

— 

0 

c 

n 

= 

0.4 

n °n 

q nq 


38 


AEDC-TR-78-1 2 


APPENDIX  A 

AERODYNAMIC  MATHEMATICAL  MODELS 

a.  Attack  Aircraft  Aerodynamic  Mathematical  Model 
(No  Rotational  Data) 

Fx  = S |cx  (n.  SH)  + ACX  (Ss.  «)] 

F'y  = S j^Cy  (a,  /3,  f5H)  -i  ACy  (SA,  a)  - ACy  (5S,  a)  ACy(5R,  a) 

+ Cy  (a)  + Cy  (a)  ii] 

1 n 2v  r 2v  J 


F ? = q^  S (c  . S,|)  - ACjr  (SA,  a)  + ACZ  (Sg,  a)  i C7  (a)  -Ip.  + C7.  (a) 

Mx  = qMSb  [Ce  (a,/3,S„)  -i-  ACj.  (SA,  a)  • ACf  (Sg,  a) 

, iCe(S„,a)  * Cf>)£  + Ct  (a)  it  + 


q 2v 


* C?A  ^’1 

lj8  2v  J 


M 


V - %.  S'  [Cm  <«■  *||>  * AC.  (SA,  «)  4C„  (8,,  a)  + ACm.  SH»  “ 

*c-,<"*»,£  *c.„  £ tC.t(A)] 

M/  = q^,  Sb  [c.  («,  /3 , <5H’>  i ACn  (SA,  a)  - ACn  (3S,  a)  - ACn  (<3p  , a) 

♦ c.r  ft  - C-P  “ ft  + \ ft  + % &] 


The  data  matrix  was  formulated  as  a function  of  the  following  variables  and  their 
associated  ranges. 

a = 0 to  90  cleg 
fi  = -90  to  90  deg 

= - 5 to  - 25  deg 


- —23  to  25  deg 
SH  = -6  to  6 deg 


5 s = 0 to  60  deg 
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b.  Fighter  Bomber  F-4  Aerodynamic  Mathematical  Model 
(With  Rotational  Data) 


FX  = S [cx  (a,  0,  aH)  + ACX  (a,  p,  5h,  nss)  + cx  (a)  J>cj 

FY  = S |Cy  + A Cy  (a,  0 , ST)  + ACy  (a,  0,  flss) 

Cy  (a)  f°_  + Cy  (a,  S.|)  —P—1 
p 2v  r 11  2vJ 

Fz  = S |cz  (a,  1 3,  8„)  + ACZ  (a , / 3,  SH,  flss)  + Cz  (a)  1°/] 


M x 


= Sb  (a  , 0)  + ACji  (a,  0,  S_x)  l-  ACj/  (a,  0, 

+ AC,  (a,  0,  il  ) + C,  (a)  .PiLb_  - C.  (a,  8„)  lA 

p 2v  T r 2v 

- c».  ii  + c,  vn 

K0  2 V lq  2vJ 

MY  = %,  Sc  |cm  (a.  0-  SH)  - ACm  (a,  /3 , «„.«„)  + Cm  (a) 

L q 2v 

* S (£) * (£)] 


M. 


= ^ Sb  [Cn  (a,  0)  + ACn  (a,  /8,  SA)  + ACn  (a,  0 , S„) 


+ AC„  (a,  0,  nj  + C„  (a)  ^ - CQ  (a,  SR)  V 


r b 


2v 


2 v 


c . ^ - cn  l°!l 

"0  2 v % 2vJ 
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The  data  matrix  was  formulated  as  a function  of  the  following  variables  and  their 
associated  ranges. 


a = 

0 to  90 

(Slcady-Statc  Data) 

55  to  90 

(Rotational  Data) 

P 

-40  to  10 

(Steady-State  Data) 

0 to  15 

(Rotational  Data) 

sn 

-2L  to  7 

5a 

-30  to  30 

5r 

-30  to  30 

^ss 

-0.3  to  0.3 
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APPENDIX  B 

ROTATIONAL  BALANCE  DATA  IMPLEMENTATION 


The  rotational  balance  data  incorporated  in  the  aerodynamic  data  matrices  of 
Appendix  A were  obtained  with  a model/balance  rotating  in  the  wind  tunnel  with  its 
rotational  vector  aligned  with  the  velocity  vector  of  the  wind  tunnel.  Because  the 
rotational  and  velocity  vectors  are  normally  misaligned  in  spinning  flight,  the  direct 
implementation  of  the  rotational-balance  data  as  a function  of  the  total  rotational  rate 
vector  (£2)  is  in  error.  The  method  used  for  implementing  the  data  resolves  a component 
(£2SS)  of  the  total  rotational  vector  along  the  aircraft  velocity  vector.  This  component 
rotational  vector  (£2SS)  may  now  be  used  in  determining  the  incremental  coefficients  to 
be  included  in  the  aerodynamic  data  matrices.  The  determination  of  the  magnitude  of  the 
steady-state  rotational  vector  (£2SS)  and  its  inclination  (6)  to  the  aircraft  total  rotational 
vector  (£2)  is  the  same  as  that  used  in  Ref.  7 and  is  outlined  below. 


Y (Body) 


V = vi  + vj  + wk 

T = pi  + qj  + rk 

^ss  = PssT  + % J + rsJ 


Total  Velocity  Vector 

Total  Rotational  Vector 

Steady-State  Rotational  Vector  along 
Velocity  Vector  v- 


|n..|  = M cosS 

-v-InlM  cos  S 
up  + vq  + wr  = |ft|  |v|  cos  8 

Substituting: 

I ft  I = np  + vq  + » 
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Knowing  the  magnitude  and  direction  of  the  steady-state  rotational  vector,  its 
components  may  be  resolved  into  body-axis  components. 

PSH  = £2SS  cos  /3  cos  a 

= fiss  sin  0 
rgs  = Qgs  cos  /3  sin  a 


These  body-axis  components  of  the  steady-state  rotation  are  now  used  for 
determining  the  magnitude  of  the  oscillatory  rate  components. 

Po  = P " P^s 

= 1 — ^ss 


The  oscillatory  components  are  used  in  conjunction  with  the  dynamic  oscillatory 
derivatives  in  the  aerodynamic  data  matrices.  Utilizing  the  rotational  and  oscillatory- 
derivatives  in  the  manner  outlined  accounts  for  oscillations  that  normally  occur  in 
addition  to  an  aircraft  steady-state  spinning  motion. 

In  the  subject  investigation,  the  rotational-balance  data  were  applied  when  the 
misalignment  angle  (5)  between  the  rotational  vector  (fi)  and  velocity  vector  (v)  was  less 
than  40  deg.  Above  40  deg,  the  rotational-balance  increments  to  the  aircraft  aerodynamic 
matrix  were  zeroed  with  an  exponential  decay  function.  Under  these  conditions,  the 
oscillatory  rates  (p0.  q0,  r0)  become  equal  to  the  total  rotational  rates  (p,  q,  r).  Likewise, 
when  the  total  rotational  vector  is  aligned  with  the  velocity  vector,  the  case  where  8 = 0, 
the  aircraft  is  in  a flat  spin  and  the  oscillatory  rate  components  are  zero. 
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b 

c 

Fx-FyjFz 

g 

lx  4y  ,h 
•xz 

Mx  ,My  fMz 
My  t 

Myg>Mzg 

m 

psq,r 

Po  >Qo  jI"o 

Pss  .qSSJrss 

P.q,r, 

q. 

s 

TxJz 

u,v,w, 


NOMENCLATURE 

Wing  span,  ft 
Wing  chord,  ft 

Forces  along  the  aircraft  body  reference  axis  in  the  X,  Y,  and  Z directions, 

lb 

Acceleration  caused  by  gravity,  ft/sec2 
Moments  of  inertia,  body  axis,  slug-ft2 
Product  of  inertia,  body  axis,  slug-ft2 

Moments  about  the  aircraft  X.  Y,  and  Z body  reference  axes,  ft-lb 

Moment  contribution  from  engine  thrust  about  the  aircraft  Y body  axis, 
ft-lb 

Engine  gyroscopic  moments  about  aircraft  Y and  Z body  axes,  ft-lb 
Aircraft  mass,  slug 

Total  roll,  pitch,  and  yaw  rates  about  aircraft  X,  Y,  and  Z body  axes, 
rad/sec,  dcg/sec 

Oscillatory  roll,  pitch,  and  yaw  rates  about  aircraft  X,  Y,  and  Z body  axes, 
rad/sec,  deg./sec 

Steady-state  roll,  pitch,  and  yaw  rates  about  aircraft  X,  Y,  and  Z body  axes, 
rad/sec,  deg/sec 

Total  roll,  pitch,  and  yaw  accelerations  about  aircraft  X,  Y,  and  Z body 
axes,  rad/sec2  deg/sec 

Dynamic  pressure,  lb,/ft2 
Wing  area,  ft2 

Aircraft  thrust  components  along  aircraft  X and  Z body  axes,  lb 
Velocity  components  along  aircraft  X,  Y,  and  Z body  axes,  ft/'sec 
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u,v,w  Acceleration  components  along  aircraft  X,  Y,  and  Z body  axes,  ft/sec2 

Free-stream  velocity,  ft/sec 

X,Y,Z  Orthogonal  body-axis  system  fixed  in  the  aircraft  with  the  X-axis  along  the 

longitudinal  centerline  of  the  body,  the  Y-axis  normal  to  the  plane  of 
symmetry,  and  the  Z-axis  in  the  plane  of  symmetry 

a Aircraft  angle  of  attack,  deg 

a Rate  of  change  of  angle  of  attack,  deg/sec 

(3  Aircraft  angle  of  sideslip,  deg 

|3  Rate  of  change  of  angle  of  sideslip,  deg/sec 

6a  Aircraft  aileron  deflection,  positive  deflection  produces  a negative  moment, 

deg 

5h  Aircraft  horizontal  stabilizer  deflection,  trailing  edge  down,  positive,  deg 

5r  Aircraft  rudder  deflection,  trailing  edge  left,  positive,  deg 

6s  Aircraft  spoiler  deflection,  deg 

9 Aircraft  angle  of  pitch,  deg 

0 Aircraft  angle  of  roll,  deg 

S2  Aircraft  total  spin  rate,  rad/sec 

nss  Aircraft  steady-state  spin  rate,  rad/sec 

I 

STATIC  AERODYNAMIC  COEFFICIENTS:  BODY  AXIS 

Cx  = Fx/q=oS  Ct  = Mx/^Sb 

CY  = Fy/q^S  Cm  = My/Oc 

Cz  = F7/qMS  Cn  = IVOb 
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DYNAMIC  STABILITY  DERIVATIVES:  BODY  AXIS 


Cy  = dCy/d  (pb/2Vso) 

p 

Cy  = dCY/dW2VJ 

r 

ct  = acj/aipb^vj 

p 

C(  = aC/(5(qc/2VJ 
q 

C„  = dCn/d(rb/2VJ 

r 

Cn^  = ^Cn/a(/3b/2VJ 


Cf  = <9C/(?(rb/2VJ 

r 

Cg.  = 3C/a(^b/2VJ 
Cm  - <?Cm/<Hpb/2VJ 

P 

= acm/d(qc/2Vj 

q 

CX  = Cx/a(qc/2VJ 
q ^ 

Cz  = dCj/diqcfflJ 

q 


C = dC  /<J(rb/2V  ) 

m m/  v oo 

r 

C - dC  /d  (ac/2V  ) 

m*  m/  v ' oo' 

a 

Cn  = acn/a(pb/2Vj 

P 

C = <9C_/<J{qc/2V  ) 

n n 1 oo 

q 
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